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ABSTRACT 

Aims. We put theoretical constraints on the presence and survival of icy grains in debris discs. Particular attention is paid to 
UV sputtering of water ice, which has so far not been studied in detail in this context. 

Methods. We present a photosputtering model based on available experimental and theoretical studies. We quantitatively 
estimate the erosion rate of icy and ice-silicate grains, under the influence of both sublimation and photosputtering, as a 
function of grain size, composition and distance from the star. The effect of erosion on the grain's location is investigated 
through numerical simulations coupling the grain size to its dynamical evolution. 

Results. Our model predicts that photodesorption efficiently destroy ice in optically thin discs, even far beyond the sublimation 
snow line. For the reference case of /3 Pictoris, we find that only <; 5 mm grains can keep their icy component for the age of 
the system in the 50-150 AU region. When taking into account the collisional reprocessing of grains, we show that the water 
ice survival on grains improves (grains down to ~ 20 /im might be partially icy). However, estimates of the amount of gas 
photosputtering would produce on such a hypothetical population of big icy grains lead to values for the O I column density 
that strongly exceed observational constraints for (3 Pic, thus ruling out the presence of a significant amount of icy grains in 
this system. Erosion rates and icy grains survival timescales are also given for a set of 11 other debris disc systems. We show 
that, with the possible exception of M stars, photosputtering cannot be neglected in calculations of icy grain lifetimes. 

Key words, stars: circumstellar matter - planetary system:formation - planetary system: protoplanetary discs - starsdndividual: 
(3 Pic 



1. Introduction 

Water ice has been unam biguously detecte d towards em- 
bedded YSOs/protostars ( Gibb et ai]|2004 . and references 
therein) and around Herbig AeBe and T Tauri stars (e.g. , 
Creech-Eakman et al.ll2002l;lMeeus et al.|[200lllGibb eta! 



20041) . For the more evolved debris disc systems, however, 



direct observational clues for the presence of icy grains 
seem far less conclusive. The obvious exception is of course 
our own Solar System where, in addition to appearing 
on satellite and planetary surfaces, dust-sized ice parti- 
cles have been detected in the comae of comets and in 
planetary rings. For extrasolar systems, the presence of 
ice appears very likely, given the fact that many debris 
discs extend to cold regions very far from their central 
stars, beyond the expected sublimation limit. However, 
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no unambiguous "smoking gun" (presumably water gun?) 
signature of their presence is available yet. Water ice is 
usually inferred through indirect arguments like best fits 
of the Spectral Energy Distribution (SED). The most con- 
clusive of these indirect detections are pro bably the ones 
around the K s tars eEridani (0. 5-1 Gyr) (|Li et al.|[2003h 
and HD 69830 (|Lisse et al.| [2006). for which the presence 
of water ice appears as a robust conclusion. 

For discs around more luminous stars, definitive con- 
clusions are more difficu l t to re ach and appear more model 
dependent. IChen et all ( 2006h analyzed spectral energy 
distributions (SEDs) around A, B and F types stars with 
detected infrared excesses and find that, when the ex- 
cesses are reproduced by a single temperature blackbody 
emission, the peak in the number of systems with a given 
grain temperature occurs at 110-120K. This was inter- 
preted by the authors as a possible indication of the icy na- 
ture of the grains, since sublimation timescales for warmer 
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grains are much shorter than the estimated system ages. 
The infrared (IR) emission detected in the HR4796 sys- 
tem (an A-ty pe star) was also reproduced by 110K black- 
body grains (jJura et al.lll998f h Absence of warmer grains 
was again interpreted as a result of sublimation of ice 
in the inner part of the disc (however, as another pos- 
sible ex planation, a sweeper com panion was also men- 
tioned). iGolimowski et al.1 (|2006t ) invoked ice to explain 
the observed color change beyond 120 AU in the archety- 
pal (5 Pic disc (again around an A-type star). Their idea 
is that beyond the sublimation zone the average grain size 
is bigger, since the grains can preserve icy mantles. 

However, even if the observational evidence is not 
100% conclusive, there are numerous theoretical argu- 
ments supporting the possible presence of ice in debris 
discs. As already mentioned, the most obvious one is that 
debris discs often extend to regions beyond the sublima- 
tion snowline, i.e. thought to be cold enough for water va- 
por to condense on dust grains, the exact location of the 
sublimation zone depending on the thermal balance at the 
grain surfaces. Furthermore, since ices have been detected 
around YSOs, it should be logical to expect it to still be 
there at later times, provided no new ice removing m echa- 
nism appeared along the way. Li fc Greenberd ( 1998a ) also 
argued that if grains are produced from cometary evap- 
oration then, since comet dust may consist of primitive 
interstellar dust on which the presence of an icy mantle 
is a well established fact, we should expect evaporated 
volatiles to recondense in the outer regions of debris discs 
(beyond ~100AU for the (3 Pic case considered by these 
authors), or the ice mantles to be preserved for dust par- 
ticles produced by shattering or collisions. 

In this paper, we revisit the theoretical constraints on 
the presence of ice in debris discs. Our main point is that 
sublimation is not the only process which leads to ice de- 
struction. We focus here on UV sputtering (a.k.a. photo- 
sputtering or photodesorption) of circumstellar ice parti- 
cles, a phenomenon well known in the context of interstel- 
lar grains. In the solar system, the ice destruction rate due 
to UV radiation was es timated to be ~ OicmMyr" 1 at 



the distance of Saturn ( Harrison! 1967t Carlson 1980l) . In 
his early pioneering work, Artvmowiczl (jl994T 19961 ) pre- 
dicted this rate to be orders of magnitude higher around 
A-type stars, which have a much stronger UV flux. In our 
paper, we provide a detailed model of the UV sputtering 
mechanism (Sect. I2.2p . It extends previous calculations 
(which were not described in detail by Artymowicz), to 
include new laboratory and theoretical results that have 
become available in the last decade, as well a range of dif- 
ferent host stars. In addressing the issue of ice survival in 
circumstellar discs, we also include sublimation fSect. l2~T|) 
in our model. In Sect. 12.51 we consider how grain dynam- 
ics is affected by erosion of the grain surface. Section 13.11 
deals with application of our study to the (3 Pic system. 
Sputtering rates for 11 other systems with detected debris 
discs can be found in Sect. 13.21 Summary and Conclusions 
are presented in Sect. [4] 
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Fig. 1. Absorption coefficients for water ice grains of dif- 
ferent sizes. (Mie theory code of Artymowicz 1988.) 



2. Model description 

In order to address the issue of ice survival in debris discs 
around stars we consider, in that order: (i) sublimation, 
(ii) photosputtering, and (iii) related grain dynamics. We 
compare the erosion timescales with the system age as 
well as the collisional lifetime of grains, the latter being 
importantly a source of newly generated dust grains. 

2.1. Sublimation 

Sublimation is a well known process. The mass loss rate 
<& is proportional to the difference between the phase- 
equilibrium vapor pressure p sa ± and the actual water vapor 
pressure p above the icy surface, i.e. <& oc (p sa t — p)u In 
our calculations, we neglect the possible presence of water 
molecules in the ambient gas, assuming p = (an assump- 
tion which is justified in a debris disc). In this case, the 
mass sublimation rate can be written as (|Lamvlll974r ) 



$ = 4.08 x 10~ 2 



Psat 

1 tort 



Zi- 



HC 



1/2 



g cm 2 s 1 , 



(1) 



where p sa t is the phase-equilibrium vapor pressure at the 
surface temperature T, and fi is the atomic weight of 
the water molecule. There are different available prescrip- 
tions, usually valid over different temperature ranges, for 
calculating the saturated vapor pressure oyer water ice 
surfaces (e.g., Marti fc Mauersbergerl Il993l Buck 1981 
ISmithsonian Meterological Tabled Il984h . For the T 



> 



170K range, we adopt the formula from lFanale fc Salvaill 
(1984) 



Psat = 2.67 x 10 10 exp 



6141. 667K 
T 



torr, 



(2) 



Condensation occurs when p > p sat . 
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For T < 170 K, we follow the experimental ly derived re- 
sults of Mauersberger fc Krankowskv ( 20031 ) 



p sat = 5.69 x 10 1J exp 



7043.51 K 
T 



torr, 



(3) 



which converges with Eq. [2] at 170 K. The erosion rate of 
a spherical particle due to sublimation is then given by 



^sub 



p 



Airps 2 



(4) 



where s is the grain radius, m su b is the rate at which a 
grain loses mass, p is the density of the grain material, i] 
is the covering factor (fraction of the surface covered by 
sublimating material). The sublimation lifetime of a grain 
then reads 



^sub — 



so 

^sub 



(5) 



where sq is the initial radius of the grain. 

Sublimation also affects the energy balance of a grain 
by removing heat at the rate 



E suh = 4tts 2 77 L s *(T), 



(G) 



where L a is the latent heat of sublimation. Strictly speak- 
ing, L s is a function of the grain temperature, how- 
ever the dependence is weak (2.6 x 1 10 erg g" 1 at OK, 
2.85 x 10 10 erg g" 1 at 90- 300 K, see lLienlTl990l ) and is 
not important for our calculations. We therefore adopt a 
value of 2.85 x 10 10 erg g _1 through out the computation. 
At T ^ 150 K, E suh is negligible compared with the energy 
loss rate due to thermal radiation (see Sect. [3]). 

Sublimation acts in the inner (warmer) part of a disc 
and destroys icy grains, or at least de-ices them, very effi- 
ciently. The radial dependence of the speed of evaporation 
process is very rapid, owing to $ varying exponentially 
with 1/T which, in turn, depends on the distance from 
the star R. This allows us to define the snowline location 
based on virtually any answer to the question of what 
constitutes a "fast enough" evaporation rate: widely dif- 
ferent definitions of that rate all result in nearly the same 
locations of the snow line (see Sect. I3.1TT|) . 

2.2. Photosputtering 

Even outside the thermal evaporation zone, grains could 
be subject to destruction through a mechanism that does 
not have nearly as sharp a spatial boundary. Individual 
UV photons absorbed by a grain do not only dissoci- 
ate water molecules but can cause molecule desorption 
from the surface of a grain. This process is known as pho- 
tosputtering or photodesorption. Around the UV-bright 
stars the effect of photosputtering can be rather strong 
(| Artvmowica 119961 ) and the position of the actual snow- 
line is no longer determined by sublimation alone. This 
is particularly true for debris discs, which are optically 
thin, so that the outer parts are not shielded from the UV 
radiation. 
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Fig. 2. Stellar fluxes, normalized to 1 AU distance from 
the star. The vertical solid line corresponds to 5.1 eV (the 
energy of the O-H bond), the vertical dashed line marks 
7 eV (the energy at which the water ice absorption coeffi- 
cients drop, Fig.[T]). 



The mechanism of photodesorp tion at the atom i c leve l 



Pt 

has recently been investigated bv I Andersson et aL ( 20061 ) 



with their theoretical modeling of the interaction be- 
tween 8-9.5eV photons and icy surfaces. First, a water 
molecule is dissociated by a UV photon. After that, in 
some cases the OH~ fragment can be directly desorbed or 
the molecule recombines and may leave the surface as a 
whole or it can pass its momentum to one of the neigh- 
boring molecules, which can leave the surface if the mo- 
mentum is in the right direction. These authors estimated 
the desorption probability, called (photosputtering) yield 
Y, and found it to lie between 4 x 10~ 4 for amorphous 
ice and 2 x 10~ 3 for crystalline ice. Y ield estimates were 
also obtained bv lWestlev et~aH (|l995h . but this time from 
laboratory experiments, studying UV sputtering of water 
ice by Ly Q photons (10.2 eV). This study found that the 
yield is a weak function of the grain surface temperature 
and of the irradiation dose. At any given fluence (flux x 
exposure time), Y is slightly higher for warmer grains. For 
a constant grain temperature, the yield increases rapidly 
with irradiation dose and reaches a p lateau value for dose s 
higher than ~3x 10 18 photons/cm 2 ( Westlev et al.lll995h . 
This plateau value is on the order of few x 10 -3 at T=35- 
100 K. The obtained temperature dependence is relatively 
weak, the extreme values being Y = 3 x 10~ 3 for T < 35 K 
and Y = 7 x 10~ 3 for T > 100 K. 

Direct compa r ison between the experimental results 
bv IWestlev et al. I dl995h and those derived from theoreti- 
cal calculations by Andersson et al. (2006) is not straight- 
forward. On the one hand, the latter estimates were de- 
rived for a fixe d and very low temp erature of 10 K. On 
the other hand, Westlev et al. ( 19951 ) data were obtained 
through laser beam irradiation, which might have induced 
local point heating and thus adding some sublimation to 
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the pure sputtering effect (N. J. Mason, private commu- 
nication) . It is thus difficult to precisely derive from these 
studies the exact dependence of the yield on the energy of 
an absorbed photon and the grain temperature. However, 
it is worth noticing that both these independently de- 
rived ranges for Y are comfortably of the same order of 
magnitude. We will thus make the simplified assumption 
that Y does not depend on temperature (which is rea- 
sonable, given a we ak temperature dependence found by 
Westlev et ail Il995l) and is furthermore independent of 
the absorbed photon wavelength The constant Y value 
we assu me is 10~ 3 , intermediat e between the lowest esti- 
m ate of lAndersson et all (|2006h and the highest estimate 



of lWestlev et a 1.1 (|l995l) 



A crucial issue for the present problem is which en- 
ergy range should be taken into account for UV sputter- 
ing calculations. There is indeed no reason why it should 
be restricted to the 8-1 leV range explore d in the 2 afore- 



mentioned studies. Dominik et al. ( 2005h . in their study 



of an origin of gas-phase water in the surface layer of pro- 
toplanetary discs, considered photons with energies down 
to 6 eV and took a constant sputtering yield per incident 
photon. However, since desorption is triggered by the ini- 
tial dissociation of a water molecule (even if the molecule 
can recombine afterwards), we adopt here a slightly lower 
threshold value of 5.1 eV, which correspond to the disso- 
ciation energy of the O-H bond ( Ruscic et al.ll200 2)FI An 
important point is that only absorbed photos can play 
rol e in the photosp u tterin g pr ocess. Photons consid ered 
by lAndersson et all (|2006h and IWestlev et all (|l995l ) be- 
longed to the 8-11 eV energy range. For such photons the 
absorption coefficient by water ice is ~ 1 (Fig.QJ, thus the 
number of absorbed photons was very close to the num- 
ber of incident photons. For this reason, in this energy 
range it does not matter whether one considers the num- 
ber of insident or absorbed photons. However, water ice 
does not efficiently absorb UV photons in the E p h ^ 8 eV 
range, and the difference between incident and absorbed 
fluxes is significant. Although for photons in the 5.1- 7 eV 
energy range Q' abB <C 1, the effect of such "low energy" 
photons cannot be neglected, because many stars radiate 
significantly more in the 5-7 eV range than above 8 eV 
(see the /3 Pic case in Fig. (2). However, the excact value 
of the lower _E p h boundary is not critical, provided it is 
below ~ 6.5 eV. 

For the upper boundary of the ene rgy of sputtering 
photons, we follow Dominik et a l. (2005) and take 13.6 eV. 
The exact value for this upper limit is not important as 
the amount of radiation beyond lleV can be neglected 
for most stars (see Fig. [5]) 

For an icy grain, the erosion rate due to UV sputtering 
is calculated then as 



(7) 



2 This energy still exceeds the surface binding energy of a 
lattice p article for water ice , which is on the order of 0.1- 



where rj is the fraction of surface covered by ice, TO(h,o) = 
3 x 10 -23 g is the mass of a water molecule and A^bs is 
the number of absorbed photons 



AUs = 



F(R,X) 
hc/X 



XX) dx, 



(8) 



in which F(X, R) is the incident stellar flux at the position 
R of the grain, A m j n = 0.091 jim (13.6 eV), A max = 0.24 ^im 
(5.1 eV), Q a bs(^) i s the absorption coefficient for pure wa- 
ter ice surfaces (Fig. [T]). 

2.3. Gas production 

Sublimation and photosputtering of ice produce gas. A 
significant fraction of gas in debris discs could be produced 
by the solids rather than be leftovers of the primordial 
dust component. This opens a possibility to deduce dust 
composition from abundances of various species in the gas 
component (if we know the parameters of the desorption 
processes) . 

In the case of icy grains, gas production rate due to 
UV sputtering can be calculated as 



Mi 



H 2 



rjAirps So 



1 AU 
R 



dn(s) dz 2nRdR 



SirpsQi] (1 AU) 2 / (^)dR, 



(9) 



where sq is the sputtering rate at 1 AU from the star, dn(s) 
is the number density of grains of size s, z is the height 
above the midplane, and t±(R) is the normal geometri- 
cal optical thickness of the disc. If a disc has a constant 
opening angle, i.e. a constant H/R ratio, where H is the 
height scale of the disc, then the above equation can be 
rewritten as 



M H2 o = SnpsoV (1 AU) 2 T|, ( — 



(10) 



where n is the geometrical optical thickness of solids in 
the disc along a radius in the midplane. 

2.4. Grain temperature 

As we have discussed, both sublimation and photosput- 
tering erosion rates depend on grain temperature. This 
dependence is weak for UV sputtering but temperature 
obviously plays a crucial role for sublimation and should 
thus be accurately computed. The grain temperature is 
calculated from the energy balance equation, in which the 
radiative heating by the central star is balanced by energy 
loss to thermal re-radiation (assuming uniform tempera- 
ture of the surface) and sublimation 



(f ) 2 /o°°Qab s ( S ,A)F,(A)dA = 

4 [Jo 00 Qabs(s, A) itB\{T) dX + n L B $(T)] 



(11) 



0.37 eV (|Diikstra et al.ll2003l ). 



where r is the distance from the star, i?* is the stellar 
radius, F* is the stellar flux at B\(T) is the Planck 
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function at equilibrium temperature T, s is the grain ra- 
dius, and Qabs is the absorption coefficient for the whole 
grain. Notice that in Eq. [TT] we use the "true" Q a b s ob- 
tained for a given grain composition whereas in Eq. [5] we 
use the absorption coefficient Q' of pure water ice. The 
reason for this is that the grain temperature depends on 
absorption by the whole grain, while the UV sputtering 
on absorption by water molecules near the the surface. 

The absorption coefficient Q a bs is a function of gr ain 
size and chemical composition fe.g. lBurns et al.lll979f ). It 
is here calculat e d usin g the Mie theory code developed 
by lArtvmowic j (1988). Effective medium theory is used 
to obtain effective dielectric constants or equivalently re- 
fraction indices, based on material component's properties 
and volume mixing ratios. Q a b s plays a crucial role in cal- 
culations of the grain temperature and, as a consequence, 
in determination of the T(R) dependence. A previously 
popular piecewise power-law ansatz (Q a bs = 1 for A < s, 
Qabs = s/X otherwise) is not suitable for calculating the 
temperature of the predominantly icy grains, because it 
gives a significant overestimate of the temperature and, 
as a consequence, much overestimated radius of the ice 
sublimation boundary. 

2.5. Icy grain erosion and dynamics 

In our model we not only estimate sputtering and subli- 
mation rates, we also compute the effect these mechanisms 
have on grain dynamics in a collisionlcss, gas-free disc. We 
follow the evolution of an icy grain taking into account the 
above two erosion modes, and the resultant change of the 
size-dependent radiation pressure coefficient (ratio of ra- 
diation pressure to gravity) . The equation of motion of a 
grain reads 



d 2 R 

dt 2 



GMr, 



■R{1-0)+F pr , 



(12) 



where m and R are the mass and position of the test 
particle, G is the gravitational constant, and M is the 
stellar mass. -Fpr is the Poynting-Robertson (PR) drag: 



PR 



-0 



GMm 



v r R 
~c~R 



(13) 



where v is the velocity vector of the grain, v r is its radial 
component, and c is the speed of light. 

The coefficient depends on the grain size, chem- 
ical composition and stellar luminosity (more precisely, 
luminosity to mass ratio L/ M). In our code, is calcu- 
lated using the routines from I Artvmowicz! ( 19881) . and dis- 
played in Fig. [3l For sizes above a certain threshold value 
(s <; 0.3 jum for the case displayed here), 0(s) is well ap- 
proximated by the geometrical optics approximation 



= 0.5 



«l/2 



(14) 
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Fig. 3. Ratio of the radiation pressure force to the gravita- 
tional force vs. grain size for different grain compositions 
around an A5V star (like Pic). Ice means amorphous wa- 
ter ice. The numbers in parentheses indicate the volume 
fraction of the corresponding elements. 



where S1/2 is the grain radius for which — 0.5. In this 
approximation, the gradual decrease in grain radius s due 
to erosion will cause slow changes in at a rate 



Sl /2 

= -O.5^s, 



(15) 



where s — s su b + s sp . Since s < 0, > 0, at least for 
grains bigger than few tenths of a micrometer (which is, 
incidentally, the wavelength of the peak power of radiation 
exerting the pressure). As can be seen from Eq. [TSJ such 
gradual changes in will increase the outward radiation 
pressure, which will tend to push the particle outwards, 
but at the same time, also increase the PR drag which 
tends to move particles inwards. For a grain on a Keplerian 
orbit of semi-major axis a, one can show that, for small 
eccentricities, the rate at which a changes reads 







20GM 



1-0 {l-0)o?c 1-0 



0» 



- t 



PR 



(16) 



where ipR = a 2 c/(20GM) is the Poynting-Robertson 
time. If > tp R , i.e., if the erosion rate is larger than 
about 2 s 2 /(si/2 ipR.), or the erosion timescale s/dots is 
shorter than (si/ 2 /2s) ipR, then PR inward migration is 
replaced by the outward erosional migration. 



3. Results 

We now present typical grain behaviors for different rep- 
resentative stellar types. For sake of clarity, we first con- 
sider a nominal case, for which we present the results in 
detail, to fully describe the mechanisms at work. For this 
reference case, we consider the archetypal Pic system. 
Results for other stellar types are presented in the follow- 
ing subsections. 
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3.1. A detailed example: Pic 

[3 Pic is an A5V main-sequence star surrounded by an op- 
tically thin debris disc, which was directly imaged for th e 
first time more than 20 years ago (|Smith fc Terrildll984l) . 
We take this system as our reference case, since it is still 
by far the most extensively studied and best known de- 
bris disc system. Countless observational and theoretical 
studies have b een dedicated to this system (e.g., the most 
recent ones: Chen et al. 2007 ; Golimowski et al. 2006: 



Thcbault & Augercau 



20061 
20051: 



Tamura et all 120061: iGalland et all 120061: iRoberge et al 
20061: Fernandez et al 



Telcsco et alj 



2005 



2004; 



Brandekcr et al 
Karmann et all 120031 : lAugereau et all l200ll ).~The disc 



has a wide radial extension (up to 1000 AU), but most 
of the detected dust is thought to reside in a relatively 
narrow region between 80 and 120 AU (see for example 
th e fit of the scattered a nd thermal light profiles derived 
bv lAugereau et al~ 2001 ). In a few studies water ice was 
considered as a possi ble grain component. As was already 
mentioned in Sect-fn TGolimowski et ah ( 20061 ) considered 
the presence of icy mantles as a possibl e explanation 
for th e observed color changes. Likewise, Tamura et al.1 
(200(1) needed an icy component to reproduce their 
high resolution near-infrared and previous polarimetry 
observ ations in terms of dus t scattering. In the same 
spirit, lAugereau et al.1 (|200lh assumed that ice might 
be present beyond the sublimation boundary. However, 
as for most debris disc systems, there is no "direct" 
unambiguous observational proof for the presence of icy 
grains. So far, the existence of ice is thus still conjectural. 

We address this issue by careful modeling of the crite- 
ria for the survival of water ice. The sublimation boundary 
alone is first considered in a "static" case when grains are 
assumed to stay at a fixed distance from the star. The 
photosputtering rate and corresponding survival times in 
the same "static" assumption are estimated afterwards. 
Finally, the coupled effect of these mechanisms is inves- 
tigated when the grain dynamical evolution is taken into 
account . In doing so, we imagine ice to be predominantly 
in H2O grain mantle, thus covering the surface (77 = 1). 
This maximizes the erosional and dynamical evolution. 

3.1.1. Sublimation boundary 

The sublimation boundary can be defined as the largest 
distance from the star at which the sublimation time of 
predominantly ic y particles is short er than any other rel- 
evant time scale ( Artvmowicd 1997 ) . We take here a very 
conservative range and consider the 2 extreme timescales 
given by the orbital period t or b, in the l-100yr range. 



and the system's age t aee , i .e. 8-20Myr ([Zuckerman et al 



200ll: |Pi Folco et ail 12004 ). The crucial parameter is the 
temperature, since t su b oc exp(— 1/T) (Eqs. Q] to [5|). The 
temperature dependence is in fact so sharp that the tem- 
perature for which t su b(T) — £ or b, i.e., ~ 140K, is rel- 
atively close to the one for which t su b(T) = < a gc, i-e. 
~ 100 K. As a consequence, we shall as a first approxima- 



le+08 
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le+06 
le+05 
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Fig. 4. Sublimation timescales for grains of different sizes 
and compositions, in the "static" assumption where grains 
are assumed to stay at a fixed radial distance. Dark blue 
lines: pure water ice grains, which are the coldest ones. 
Light brown lines: icy grains with 10 % of graphite in- 
clusions, which represent some of the "hottest" grains. 
These two choices of grain compositions are not the most 
representative from the point of view of grain chemistry. 
However, they bracket temperature values for a grain of a 
certain size at a give distance from the star. 



tion neglect the timescale issue and conservatively con- 
sider that the sublimation boundary lies somewhere in 
this 100- 140 K range. The temperature of a given grain 
depends schematically on 3 parameters: its size, chemi- 
cal composition and distance to the star. This is illus- 
trated in Fig. [5] where lOOK-isotherms, giving the outer 
limit of the sublimation boundary, are plotted for differ- 
ent grain compositions as a function of grain size. If we 
consider a reasonable intermediate compositional case of 
dirty ice grains, we see that the sublimation boundary lies 
somewhere between 5-10 AU for large mm-particles and 
20-50 AU for hotter submicron-scale grains (Fig. [6|) . For 
compact grains these values are all far inside the location 
of the color change observed bv iGolimowski et al. ( 20061 ) 
at ~ 120 AU. Only very small 0.1 fim grains of special 
chemical compositions have an outer sublimation bound- 
ary approaching 120 AU (Fig. [5]), but those grains should 
be blown out by radiation (/? > 1). Thus, in the case 
of compact grains, our results rule out sublimation as a 
simple explanation for a hypothetical compositional tran- 
sition at ^120 AU from the star. Howe ver, if grains are 
highly porous, as e.g. was introduced in iLi fc Greenbergl 
(1998b). and consist of many individual submicron grains, 
the result can be different, since such aggregates optically 
behave like invidual submicron grains. 
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Fig. 5. Isotherms (T = 100 K) as a function of grain size 
for different grain compositions. These lines provide the 
outer limit for ice sublimation boundary. This plot clearly 
shows how the sublimation boundary depends on the grain 
sizes and chemical compositions. The horizontal line is for 
the blackbody approximation. The solid black line is for 
the artificial piece-wise power low dependence ( Q a b s = 1 
for A < s, Q a bs = s/A otherwise). All other lines are the 
same as in Fig. [3] 



s=0.1 um 

1 um 

■ ■ ■ ■ 10 urn 

100 um 

1 mm 




10 
R[AU] 



1000 



Fig. 6. Temperatures of dirty ice grains (60% of water ice 
+ 40% of MgSiOs) as a function of the distance from the 
star (for the f3 Pic case). The sublimation boundary lies 
somewhere in the crosshatched region. The grain radii are 
indicated on the plot. 

3.1.2. Recondensation 

Here we justify that the recondensation rate is very small 
and can be safely neglected in a debris disc of a late type, 
nearly gas free system like, e.g., (3 Pic. 

The condensation rate can be computed with the help 
of Eq.[T]if instead of p sa t, the actual water vapor pressure 
above the ice surface p is substituted. In the ideal gas 



Table 1. Sputtering rate by UV photons at 1AU from 
(3 Pic for different grain radii. 



s (um) s sp (/um yr ) 



0.1 
0.3 
1 

3 
7 

> 10 



1.3 
1.8 

2.35 
3.23 
4.14 
4.45 



approximation, p = kn}i 2 oT g where k is the Boltzmann 
constant, nn 2 o is the water molecule volume density in 
gas phase, T g is the gas temperature. If the gas and the 
grain surface are in thermal equilibrium T g = T, the con- 
densation rate, «3> c , can be calculated as 



$ r = 1.8 x 10" 



-26 ( n H 2 o 



1 cm 



-3 



T 

Ik 



1/2 



— 2 — 1 

g cm s . 



(17) 



To estimate 4> c we need to know ?ih 2 o- The upper 
limit can be obtained from the H distribution given in 
Brandeker et al.1 (|2004t ) for the solar composition case un- 
der an oversimplified assumption that all oxygen atoms 
are locked in water vapor. The assumption is not real- 
istic, since some fraction of oxygen is present in atomic 
phase, some is locked in CO and O2, but it gives an easy- 
to-estimate upper limit. This estimate^ gives us ?t.h 2 o < 
1 cm -3 , resulting in <E> C ^ 10 25 g cm ~ 2 s^ 1 and a reconden- 
sation rate that exceeds the sublimation rate at T ^ 85 K. 
Thus recondensation dominates sublimation in the outer 
regions of the disk, but the time scale, even for the case of 
orders-of- magnitude overestimated H2O vapor densities, is 
so high (it would take about 10 12 yr to accumulate 1 /im 
layer of ice on a grain surface) that this process can be 
safely ignored in the frame of debris disk systems. 

3.1.3. UV sputtering 

Since the disc is optically thir0, UV sputtering occurs even 
in the outermost regions. Taking into account the geo- 
metrical dilution of the stellar UV radiation and adopt- 
ing our conservative and temperature-independent value 
of Y = 10~ 3 (see section [272]) . we obtain for the UV sput- 
tering rate 



-4.45 



1 AU 
R 



fim yr 



(18) 



Detailed calculations of the ionization disk structure and 
the equilibrium chemistry network, show t hat fo r the so- 
lar composition case from iBrandeker et al.l (|2004f ) tih 2 o 
10~ 6 cm" 3 and in the metal depleted case < 0.1 cm" 3 (Liseau, 
R. private communication) 
4 Neither the amount of dust nor the a mount of hydro- 



gen ( the upper limits: iV(HI)<^ 5 x 10 9 cm (IFreudling et al 



1995 ). AT(H 2 )£ 3 



(|Lecavelier des Etangs et al 



200 if ) 1 ) is enough to effectively shield the UV radiation 
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Fig. 7. Thick black lines: UV sputtering destruc- 
tion timescales for icy grains at different distances from 
(3 Pic central star (for Y(T) = 1CT 3 ). The calculations 
are done for the "static" case, i.e. grains are assumed to 
stay at the same distance from the star. Thin red lines: 
Removal timescales i rem derived from Fig. [9] when taking 
into account grain dynamics. 



The front coefficient is a weak function of grain size 
(Table [l}, since the number of absorbed photons de- 
pends on Q^bs- The corresponding sputtering destruction 
timescale i sp = so/s sp has been plotted in Fig. [7] A first 
important result is that the dependence with R, while be- 
ing significant, is much less pronounced than for the sub- 
limation case (Fig. 3]). There is no longer an abrupt tran- 
sition between the non-eroding and eroding regions, and 
it is difficult to define a "snowline" due to sputtering. A 
consequence of this result is that sputtering is unlikely to 
explain sharp observed changes such as the color transi- 
tion at 120 AU. 

Another important result is that no "dust" grain (i.e. 
in the Si 1-5 mm range) can survive in the inner R <, 
200 AU region over the system's age. Even in the outer 
regions, only relatively large grains, in the sub to mil- 
limetre range, can survive. We shall consider, as a first 
approximation, that the minimum size for a grain which 
can retain its icy component over the system's age in the 
~80-120 AU region where most of the f3 Pic dust has been 
observed is ~ 5 mm. 



3.1.4. Icy grain dynamics 

All the time scales previously discussed were given for the 
"static" case, when a grain was assumed to stay at a given 
distance from the star. This is of course a simplification 
of the behavior of a "real" particle. Indeed, as a grain 
gets eroded, be it by sublimation or UV sputtering, its 
response to radiation pressure (measured by its (3 param- 
eter) changes. As has been shown in Sect. I2.5( this has 
consequences on the grain's orbital parameters, leading to 
outward or inward migration depending if the ratio be- 
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Fig. 8. Dynamics of UV sputtered icy grains, initially 
placed at 100 AU from the star at a keplerian velocities. 



tween the size erosion timescale is shorter or longer than 
the typical Poynting- Robertson timescale £pr. 

The full dynamical effect of size erosion has been taken 
into account in a batch of simulations where test particles 
motions are obtained by integrating Eq. [T2] (Fig. [5J . Our 
results show that for all initial particle sizes considered, 
we are here always in the outward migration regime. All 
particles released at 100 AU reach 1000 AU on a timescale 
comprised between a few 10 4 yr for the smallest considered 
bound grains (4 /nm) and a few 10 5 yr for large submil- 
limetre grains. An interesting consequence of these strong 
outward migrations is that no grain ever gets fully eroded 
to s = 0. In fact, once grains enter the outbound orbit 
regime (f3 iS 0.5), the outward migration rate increases 
very sharply and gets much higher than the size erosion 
rate. This imbalance between migration and erosion rates 



A. Grigorieva et al.: Collisional avalanches in dusty discs 



9 



le+07 



le+06 



£ le+05 




S le404 fl^T' 



1 |J.m " 

■ ■ ■ 10 p.m 
-■- 20(J.m 
--"lTjo"|im 

1 mm 

— » 5 mm 



50 



100 150 
R[AU] 



200 



250 



Fig. 9. Sputtering removal timescales t rcm and collisional 
lifetimes f co u for icy grains of different sizes. All t Tem curves 
appear as black lines, all t co u as orange ones. The increase 
of £ co ii inside 80 AU is due to depletion of dust in this 
region. 



is self-amplifying: as the grain gets ejected its erosion rate 
rapidly decreases (since t sp oc R 2 ). Eventually, the grain 
gets ejected from the system before it has time to be 
fully eroded, and its size no longer evolves once the grain 
has passed beyond ~ 1000 AU and reaches an asymptotic 
value (see Fig. [5fc). Interestingly, this asymptotic value, 
which is in the 1-2 /im range, is relatively independent of 
the grain's initial size sq, the only difference being that it 
is reached faster for initially smaller particles. 

Nevertheless, the end result is the same as if the par- 
ticles had been fully eroded: they are all eventually re- 
moved from the system. Moreover, for the bigger grains, 
the timescale i rcm for this removal (which we conserva- 
tively define as the time it takes for a grain at initial 
distance distance Rq to migrate beyond 2Rq) is compa- 
rable to t sp estimated in the previous section. Indeed, for 
a grain in the submillimetre to millimetre range, most of 
the sputtering erosion takes place at a constant distance 
from the star. It is then followed by a relatively brief ejec- 
tion phase, which begins once its size has dropped below 
the (3 w 0.5 limit, so that the total time from start to 
ejection is almost equal to the erosion time t sPl to reach 
= 0.5 (si/2 = 4/im). This time t sPi/2 is in turn almost 
equal to the full erosion time t sp (the erosion time from 
4/zm to being negligible compared to t sPi/2 ). For grains 
with smaller initial sq, the situation is less simple, the time 
required to reach Si/ 2 — 4/im getting significantly shorter 
than i sp when sq is close to Si/2- Furthermore, the ejection 
timescale t cj is no longer negligible compared to t SPl and 
has to be taken into account when estimating i rom - For 
even smaller unbound grains, of course, sputtering plays 
only a very minor role. These particles are ejected before 



A set of additional runs has been carried out, where 
t Tem is estimated, for the 1 /^m to 5 mm size range, as a 
function of initial distance from the star (Fig. [9| . As ex- 
pected, differences with the "static" case increase as so 
gets smaller. For the smallest grains, i rcm is much shorter 
than t sp and is comparable to the orbital timescale at 
100 AU, i.e. - 10 3 yr- 

3.1.5. Comparison with collision timescales 

In Fig. [9l we present a summary of all our results regarding 
ice survival times in the (3 Pic system. A first conclusion 
is that in the region where most of the dust is observed 
(80 R Ss 150 AU), sputtering is always the dominant ice 
removal mechanism. We also clearly see what has been al- 
ready mentioned in the previous subsections, i.e. that for 
all particles in the dust size range, ice survival timescales 
ij ce (be it due to sublimation or sputtering erosion com- 
bined to dynamical ejection) are always shorter than the 
system's age. We confirm the result derived in the static 
case, i.e. that the smallest icy grains which can survive 
in the whole > 80 AU region over ~ 10 7 yr have a size 
Ssurv — 5 mm. 

However, the system's age is probably not the only 
relevant timescale, t- lcc has to be compared to. Indeed, 
it is well established that no debris disc grains are pri- 
mordial, but that they should be continuously resupplied 
by collisions among b igger and undetectable objects (e.g., 
IThebault et alJl2003L for a detailed discussion on the sub- 
ject). In this small and (partially) icy grain could 
be produced from a bigger particle for which t lcc > t agc 
and possibly be destroyed by another collision before sput- 
tering is able to fully erode it[f| 

A first issue is here to see if collisions cannot by them- 
selves evaporate all ices from a grain's surface. From the 
estimates given in Table I of iTielens et al.l <| 1994T ) , the 
threshold impact velocity for achieving vaporization of wa- 
ter ice is i>th — 6.5 km s -1 . This value is very high: it ex- 
ceeds the orbital velocity beyond ~ 25 AU, thus making 
the occurrence of collisions at v > vth very unlikely (ex- 
cept possibly for very small unbound ((3 > 0.5 grains). 

The rate at which dust grains are produced and de- 
stroyed by collisions should in principle be estimated 
with proper numerical simulations which exceed by far 



the scope of this paper (see e.g. IT hebault et al.1 12003 



Krivov et all l2006t IThebault fc Augereaull2007i) . We will 



here consider the collisiona l lifetime ^ semi-empirical for - 
mula numerically derived bv IThebault &: Augereaul ( 2007 ) 



any significant erosion takes place, so that t r , 



5 We here consider "erosion" in the extended sense of erosion 
to the limit followed by dynamical ejection. 

6 The collisional lifetime being defined as the time it takes 
for a given particle of size s to lose 100% of its initial mass, be 
it by accumulation of numerous partially eroding (cratering) 
impacts, or by one violent fragmenting (i.e. total shattering) 
event. If the system is in collisional equilibrium, t co n is also 
equal to the production lifetime, i.e. the time it takes to pro- 
duce (by collisions among bigger bodies) the total number of 
size s particles present in the system. 
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for extended debris discs, valid in the s SS 10 6 s 1 / 2 size 
range: 



^coll — icollO 



2.7 



for s < S2 



with si = \.1s\j2 and si — 100si/2, and 

0.3 

tcoll = icollO 



«2 



for s > S2 



(19) 



(20) 



where t co m = (rfi) -1 is the usual, and often misleading 
reference collision timescale directly derived from the ver- 
tical optical depth r (fl being the angular orbital velocity). 
For the radial di s tribut ion of the optical depth, we follow 
I Grigorieva et alJ (120071) and assume the best fit obtained 
bv lAugereau et all (|20o"l . 

As can be seen from Fig. HI the situation is a com- 
plex one. To better understand the role of collisions, let 
us here make the simplifying but reasonable assumption 
that a particle of a given size s is mostly produced from 
collisions involving pa rent bodies which ar e typically of 
size Sparcnt — 2 — 5s ( Thebault et al. 2003). Let us now 
consider as a reliable starting point the smallest parti- 
cles that are able to retain ices over the system's age re- 
gardless of the additional effect of collisions, i.e. the ones 
with ii CC > i ago . As already discussed, this corresponds 
to objects with size s surv — 5 mm. Fragments produced 
from these parent bodies, which are predominantly in the 
~ 1 mm range, should thus initially contain ices. Since for 
these 1 mm objects t co \\ < t lcc almost everywhere in the 
disc (beyond > 40 AU), we can assume that most of these 
grains will remain icy over their collisionally-imposed life- 
time. If we carry on with this iterative procedure, we see 
that we can work our way down to s ~ 20 /jm before t co ii 
begins to get comparable to t- lcc in some regions of the sys- 
tem. For all sizes below this limit, ice survival times are 
smaller than collisional lifetimes, implying it is unlikely 
that collisionally produced fragments can remain icy. As 
a consequence, 20 (im may be considered as the approx- 
imate lower boundary for grains which can have an icy 
component in the (3 Pic disc, at least in the 50-150 AU 
region where most of the dust has been detected. 

Of course these results should be taken with great 
care, since this very complex issue should in principle be 
addressed using full scale collision evolution simulations. 
Such detailed study will be the purpose of a forthcoming 
paper. 



Table 2. UV sputtering rates and grain survival times, 
numerically estimated for a sample of 11 debris disc sys- 
tems. All values are given at a reference distance of 100 AU 
from the central star, which is a characteristic scale of the 
density maximum in debris discs. The yield per incident 
photon is taken to be Y = 10~ 3 x Q^bs- Systems marked 
with * have their sublimation boundaries outside 100 AU. 
Thus the values given in the table do not reflect the grain 
lifetimes at 100 AU, however they can be scaled to farther 
distances. 



Name 


Type 


Ssp(pniyr 




t(yr) 








20 urn 


1 mm 


HD66591* 


B4V 


7 


(2)* 


(140)* 


HD142165* 


B5V 


0.9 


(20)* 


(10 3 )* 


Vega* 


A0V 


0.09 


(250)* 


(10 4 )* 


HR4796A 


A0V 


0.04 


600 


2xl0 4 


Fomalhaut 


A3V 


3xl0~ 3 


7xl0 3 


3xl0 5 


/3 Pic 


A5V 


4.5xl0~ 4 


5.5xlC 


4 2.3xl0 6 


HD 40136 


F1V 


5xl0~ 5 


5xl0 5 


2xl0 7 


HD 1581 


GOV 


10~ 6 


2xl0 7 


10 9 


HD 20807 


G2V 


10" 6 


2xl0 7 


10 9 


e Eridani 


K2V 


2xl0" 7 


10 8 


5xl0 9 


AU Mic 


MlVe 


4xl0 -7 


5xl0 7 


2.5xl0 9 



Using the assumptions of Se ct. I3.1.3 | together with the 



(2001). 



wo 



dust distribution derived by lAugereau et al 
estimate that, should most of the grains be icy (r/ = 
1), the corresponding production rate of gas would be 
Mn 2 o ~ lO^Mgyr -1 , and the column density produc- 
tion of oxygen No ~ 10 14 cm~ 2 yr -1 . Since oxygen is 
expected to be neutral, it is a lso dynamically inert in 
the disc ( Fernandez et al. l2006h . This means that oxy- 
gen released is not easily removed, and will only slowly 
redistribute itself in the disc. By comparing to the ob- 
served O I column de nsity iVoi = (5.5 ± 2.5) x 10 15 cm~ 2 
(jRoberge et al 1 l2006l) . we find that all the observed oxygen 
would be produced from photosputtering on an extremely 
short timescale, Ss 100 yr. Together with the fact that dust 
grains unexposed by radiation are continuously produced 
in debris from collisions, this implies that the parent bod- 
ies must be very poor in water ice. Otherwise, the effi- 
ciency of the photosputtering process would quickly refill 
the system with unrealistically large amounts of gaseous 
oxygen. 



3.1.6. Gas production rate 

For the specific case of (3 Pic, these studies of photosput- 
tering and its consequences can be taken one step further 
than estimating grain sizes which could be icy. Indeed, 
as mentioned in Sect. 12.31 sublimation and photosputter- 
ing do not simply "remove" ice from the system, they 
transform it into gas. Water molecules desorbed from the 
surface are quickly dissociated i n gas phase by the ste l- 
lar UV radiation into H and O ( Li &: Greenberd 1998a ) . 



3.2. Other systems 

Performing similarly extensive studies for all debris disc 
systems is out of the scope of the present paper. We will 
here present, for a selection of 11 representative systems, 
estimates of s sp at a reference distance of 100 AU, as de- 
rived from Eq|7] Furthermore, for each system we present 
the numerically derived values for t lcm , for a "big" 1 mm 
and a "small" 20 /im grain. As can be seen from Table [2j 
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s sp varies by 7 orders of magnitude depending on the stel- 
lar typeQ- 

A can be clearly seen, for all considered stars, the sput- 
tering removal timescale of a 20 /im grain is always shorter 
than the age of the system. The only exception is the M 
star AU Mic, where the erosion timescale is ~ 50 Myr, i.e. 
at least 4 times the age of the system (12 Myr). For big- 
ger, 1 mm grains, the situation is more favorable, at least 
around the less luminous G, K and M stars, where the ero- 
sion timescale of a 1 mm grain exceeds 1 Gyr. However, for 
F stars and upwards, millimetre-sized icy particles cannot 
survive for more than l-2xl0 7 yr, which is less than the 
minimum age of a typical debris disc. 

The above numerical estimates are presented to give 
a first estimate of photosputtering effects as a function 
of stellar type and should be interpreted carefully. They 
are given at a reference distance of 100 AU, whereas dust 
might be observed significantly closer or further away from 
the star. This limitation is, however, not too crucial since, 
as a first approximation, the erosion rates and timescales 
we computed can be rescaled with distance following a 
quadratic dependence of irradiation on R. Another issue 
that has not been addressed here is the complex inter- 
play between sputtering and collisions. Since for most of 
the presented systems, dust size and spatial distribution 
is not as well constrained as for ft Pic, such an exploration 
would be somewhat tentative. It should, however, be kept 
in mind that intense collisional activity could affect grain 
survivability. Furthermore, the microscopic structure of 
the evaporating grain will certainly undergo change in the 
course of radiative de-icing. The resultant sputtered grain 
may be much more porous than in the initial phase or, on 
the contrary, develop a microcometary structure in which 
dust is mixed with ice in the large core, and with an altered 
or ganic "crust" near the surface (a possibility mentioned 



Artvmowiczl 1997). How the mechanical properties of 



a grain evolve with time in that case and how the phe- 
nomenon of dust avalanches is enhanced by icy nature of 
grains are new and potentially important questions. 

Our preliminary conclusion is that, with the possible 
exception of M stars, there is no system where photo- 
sputtering effect can be fully neglected. Of course, each of 
these systems has to be more carefully investigated, and 
our sputtering erosion estimates utilized in more advanced 
global dust evolution models. 



4. Summary and Conclusion 

We have reconsidered the issue of icy grain presence and 
survival in debris discs by quantitatively exploring the role 
of UV sputtering in addition to the well known sublima- 
tion mechanism. Taking the archetypal ft Pic debris disc 
as a reference system, we show that sputtering is able to 



UV fluxes were obtained from the MAST archive 



http:/ /archive.stsci.edu/index.html| The systems with 
debris discs were selected from the Debris Disc Database 



significantly erode icy grains far beyond the sublimation- 
imposed snow line. Using a dynamical model coupling 
grain erosion and orbital evolution, we find that, in the 50- 
150 AU region where the bulk of ft Pic's grain population 
is thought to reside, all grains smaller than a threshold 
size Scro ~ 5 mm (basically corresponding to all particles 
detectable by observations) cannot retain their icy com- 
ponent over the age of the system. Taking into account 
possible collisional activity slightly improves the situation. 
There is indeed a non negligible size range below s oro for 
which collision timescales are shorter than sputtering ero- 
sion times, so that steady collisional reprocessing starting 
at grains just above s eTO (in the icy "reservoir" size range) 
could maintain a significant amount of smaller icy grains. 
However, when checking this icy grains hypothesis against 
the observational constraints on gas abundances, we find 
that the rate at which O I should be steadily produced by 
photosputtering would very quickly lead to column den- 
sities exceeding by several orders of magnitudes observa- 
tional upper limits. This seems like a relatively robust ar- 
gument for ruling out the assumption that, for the ft Pic 
system, the observed grain population is predominantly 
icy. 

We have explored the amplitude of the sputtering 
mechanism for other debris discs systems, deriving sput- 
tering rates and survival timescales, for different grain 
sizes, as a function of spectral type We have shown that 
for all stars other than M dwarfs, no icy grain in the 20 (im 
range can survive for more than 10 7 yr, while millimetre- 
sized ice grains can survive in G stars and below. These 
derivations are of course first order approximations. The 
obtained erosion rates have to be included in more sophis- 
ticated dust evolution models, taking for instance into ac- 
count the coupled effect of collisions and sputtering. This 
will be the purpose of a forthcoming study. 
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